arXiv:hep-lat/9701005vl 9 Jan 1997 


Preprint HLRZ 84/96 

Gauge-ball spectrum 
of the four-dimensional pure U(l) gauge 

theory 


J. Cox, W. Franzki, J. Jersak 

Institut fur Theoretische Physik E, RWTH Aachen, Germany 


C. B. Lang 

Institut fiir Theoretische Physik, Karl-Franzens-Universitdt Graz, Austria 

T. Neuhaus 

Niels Bohr Institute, Univ. of Gopenhagen, Denmark 

P. W. Stephenson 

DESY, Zeuthen, Germany 


Abstract 

We investigate the continuum limit of the gauge-ball spectrum in the four-dimensional 
pure U(l) lattice gauge theory. In the confinement phase we identify various states 
scaling with the correlation length exponent v ~ 0.35. The square root of the string 
tension also scales with this exponent, which agrees with the non-Gaussian fixed 
point exponent recently found in the finite size studies of this theory. Possible sce¬ 
narios for constructing a non-Gaussian continuum theory with the observed gauge- 
ball spectrum are discussed. The state, however, scales with a Gaussian value 
n ~ 0.5. This suggests the existence of a second, Gaussian continuum limit in the 
confinement phase and also the presence of a light or possibly massless scalar in the 
non-Gaussian continuum theory. In the Goulomb phase we find evidence for a few 
gauge-balls, being resonances in multi-photon channels; they seem to approach the 
continuum limit with as yet unknown critical exponents. The maximal value of the 
renormalized coupling in this phase is determined and its universality confirmed. 
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1 Introduction 


Since studies of lattice gauge theories began it has been conjectured that the 
pure U(l) gauge theory in four dimensions (4D) might possess two contin¬ 
uum limits at the phase transition between the conhnement and the Coulomb 
phases [1,2], provided the transition is of 2"^'^ order. The continuum limit de¬ 
pends on the phase in which it is approached. In the conhnement phase various 
gauge-balls (GB) with hnite masses, analogous to the glue-balls of pure QCD, 
as well as a hnite string tension a, are expected. In the Coulomb phase with a 
massless photon dilute magnetic monopoles might exist, whereas in the con¬ 
hnement phase they condense. These phenomena were clearly observed for 
a hnite lattice cutoh 1/a long ago [3,4]. However, the continuum limit has 
remained elusive for the subsequent 13 years. 

The main reason for this delay has been the two-state signal at the phase 
transition observed on hnite lattices with periodic boundary conditions [5,6]. 
However, it has recently been demonstrated that if, instead of such toroidal 
lattices, those with trivial homotopy group are used, the two-state signal dis¬ 
appears for the Wilson action [7], extended Wilson action (dehned below) at 
7 < 0 [8,9], and the Villain action [10]. Furthermore, in high precision simu¬ 
lations on specially constructed homogeneous spherical lattices, the hnite size 
scaling (FSS) behaviour turned out to be well described by the leading term 
of the critical behaviour with a non-Gaussian value of the correlation length 
critical exponent, v = 0.365(8) [8-10]. 

Within the intrinsic uncertainties of numerical evidence this demonstrates that 
the phase transition is of second order in some parameter region and the cor¬ 
responding continuum limit is governed by a non-Gaussian hxed point. The 
expected continuum limit is thus presumably a nontrivial quantum held theory 
which is not asymptotically free. Until now no such nonperturbatively dehned 
continuum quantum held theory has been established in four dimensions. Fur¬ 
ther investigation of this theory is therefore of high theoretical interest. 

The use of spherical lattices has been crucial for the study of the model by 
FSS methods, because this method combines scaling and hnite size phenom¬ 
ena. Thus the latter should not be distorted by topological ehects. However, 
once the critical properties are established, there is no urgent need to use 
such complicated lattices for investigations of the spectrum of the model. The 
topology of hnite lattices of volume V ought to be irrelevant for the thermo¬ 
dynamic limit, V oo, if that is taken before the continuum limit, when the 
lattice constant a —0. The restrictions imposed by the two-state signal on 
toroidal lattices can be avoided by using large V and choosing the coupling 
parameters outside the metastability region on these lattices. In fact, fairly 
precise values of u = 0.28 — 0.42 have been obtained by this approach in the 
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past. 


In this work we therefore return to the 4D toroidal lattices and adapt the latest 
methods of glue-ball measurements in QCD on such lattices to the U(l) gauge 
group. Since the most recent studies have indicated universality of the scaling 
behaviour [ 8 - 10 ], we have chosen an action reducing the two-state signal while 
keeping the autocorrelation time reasonably short: the extended Wilson action 
at the double-charge coupling 7 = — 0 . 2 . 

We have investigated in both phases on large lattices the correlation functions 
in 20 channels with the zero momentum and in some channels with the smallest 
nonzero momentum as well. We have also determined the static potential from 
the same runs. 

In the conhnement phase, we have obtained masses rrij (inverse correlation 
lengths) in most channels. Within the considered parameter range the hnite 
size effects are negligible. Also the string tension a has been determined from 
the values of the potential at the largest separations. We investigate the scaling 
behaviour of these observables, i.e. their vanishing in lattice units, when (3, 
the standard Wilson coupling, approaches its critical value jSc- 

Our main result for the conhnement phase in the vicinity of the critical point 
Pc is the evidence for two groups of GB masses with distinctly different scaling 
behaviour when the phase transition is approached, r = \P — Pd 0. Most 
of the GB masses, and approximately also ^/a, scale proportional to the 
value of the correlation length exponent being non-Gaussian, 

z/„g = 0.35(3). (1) 


This agrees with the exponent u = 0.365(8) found in FSS studies on spherical 
lattices [ 8 - 10 ]. 

However, the mass of the = O^"*" GB scales as with 

"* = 0.49(7). (2) 

This correlation length exponent agrees with the Gaussian value u = 1/2. The 
critical point, common for both groups, is at 

Pc ~ 1.1607(3). (3) 


The errors indicated in the above three equations include both statistical and 
systematic uncertainties of our results. Since the latter are not independent. 
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we also give the result for the ratio of both exponents, 


^ = 0.71(8). (4) 


These hndings suggest that in the conhnement phase two different continuum 
limits are possible. One is non-Gaussian, in which all states scaling as 
are obtained. Numerous GB states with a somewhat degenerate spectrum, as 
well as the hnite and nonvanishing string tension should be expected. Our 
recent [8-10] and present results suggest that the observables scaling as 
correspond to a nontrivial continuum theory which is not asymptotically free. 
This theory contains a light, possibly massless scalar. 

The other continuum limit in the conhnement phase, based on scaling like 
is presumably Gaussian. Of all the reliably measured GB masses only that of 
the O"*"^ state has this behaviour. 

The Gaussian exponent (2) has not been detected in the hnite size scaling 
approach [8-10]. However, the existence of two eigenvalues of the linearized 
renormalization group transformation matrix, consistent with the values (1) 
and (2) of the critical exponents, was observed back in the eighties in Monte 
Garlo renormalization group studies of the model [11,12]. These exponents can 
now be interpreted in terms of the scaling of physical observables. 

In the Goulomb phase we hnd evidence that in the vicinity of the phase tran¬ 
sition, apart from the massless vector (“photon”) state, also massive states 
are present. They appear to scale at the phase transition. We identify them as 
resonances in two-photon and three-photon channels, respectively, with sig- 
nihcant hnite size effects. As these are not yet fully under theoretical control, 
the critical exponents could not be determined. Vanishing of the resonance 
masses when [3c is approached would imply the existence of a continuum limit 
with some dimensionful parameter in this phase. 

In the Goulomb phase we also obtain from the potential the renormalized 
coupling Or = The value of this coupling at the critical point /5c, = 

0.19(1) is consistent with that in other formulations of the pure compact U(l) 
gauge theory. This is further support for the conjectured universality of the 
renormalized coupling of this theory in the Goulomb phase as one approaches 
the critical point [13-15]. 

The results are presented as follows: In the next section we dehne the model 
and describe its critical properties. In sec. 3 we discuss the results for the static 
potential in both phases. Technical details about the GB measurements are 
collected in sec. 4. In sec. 5 we then present the results for the GB masses in 
the conhnement phase. Possible interpretations of their scaling behaviour are 
discussed in sec. 6. Sec. 7 is devoted to the results of GB measurements in the 
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Coulomb phase. The conclusions and a long list of interesting open questions 
are given in sec. 8 . The appendix contains tables with more detailed data. 


2 The pure U(l) lattice gauge theory 

2.1 Action and phase diagram 


We consider pure compact U(l) lattice gauge theory with an extended Wilson 
action, 


5 = — ^ [/5cos(0p) + 7 cos(20p)]. (5) 

p 

Here, 0p G [0,27r) is the plaquette angle, i.e. the argument of the product 
of U(l) link variables around a plaquette P, and f3 and 7 are the single and 
the double charge representation couplings, respectively. Taking 0p = a'^gF^u, 
where a is the lattice spacing, and P+A'y = l/g"^, one obtains for weak coupling 
g the usual continuum action S = d^^xF^^. 

This lattice gauge theory has a line of phase transitions between the strong 
coupling conhnement phase and the weak coupling Coulomb phase. Its position 
for the Wilson action (7 = 0) is {3c — 1.011 [ 6 ], whereas for 7 = —0.2 it is 
given in (3). For 7 > +0.2 the transition is clearly of 1®* order, weakening with 
decreasing 7 [16,6]. 

The recent studies on spherical lattices strongly suggest that the order changes 
at 7 = 7 o — 0 , 7 o probably being slightly positive, and is of 2 °'^ order for 
7 < 7 o [8,9]. With decreasing 7 the 2^^^ order transition further weakens in 
the sense that the specihc heat peak decreases for hxed lattice size and the 
autocorrelation time increases [9]. The scaling behaviour of bulk quantities is 
universal at least in the range —0.5 < 7 < 0. 

On toroidal lattices the disturbing two-state signal weakens with decreasing 
7 , but is present at least until 7 = —0.5 [ 6 ]. For even smaller 7 the large 
autocorrelation time makes simulations prohibitively expensive. Thus on the 
toroidal lattices the two-state signal on the critical line cannot be avoided. 


2.2 Some earlier studies 


It has been suggested [ 6 ] that the point 70 is a tricritical point (TCP). However, 
we would like to point out that the TCP’s are not dehned merely as points 
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on phase transition lines where the change between the 1 ®* and 2 "^^ order 
takes place. More is required: in a parameter space enlarged by an external 
coupling to the order parameter three critical lines should emerge from the 
TCP [17,6]. This results in TCP’s being associated with tricritical exponents 
different from the critical ones on the emerging critical lines. Up to now there 
has been no evidence that the point 70 in the pure compact U(l) gauge theory 
is really a TCP, because the other critical lines are not known. The problem 
is in hnding the suitably enlarged coupling space while only nonlocal order 
parameters distinguishing between the conhnement and Coulomb phases, like 
the string tension or the photon mass, exist. This problem may be specihc for 
gauge theories. It makes an interpretation of the two observed scaling laws in 
terms of critical and tricritical scaling behaviour uncertain. 

Various investigations of the scaling behaviour of the model by analytic means 
[18], the hnite size scaling analysis [16,19] and the Monte Carlo renormalization 
group (MCRG) method [11,12,20] found the correlation length exponent v to 
be in the range v ~ 0.28 — 0.42. In some MCRG studies the existence of a 
second value of this exponent, consistent with v = 0.5, was observed [11,12] 
suggesting that one of these exponents is tricritical and the other critical. In 
ref. [11] the tricritical value was suggested to be z/ = 0.5, whereas in refs. [12] 
the assignment was opposite, the tricritical exponent being the non-Gaussian 
one. Recently, evidence has been provided, that in the coupling space enlarged 
by monopole coupling, a second order phase transition with non-Gaussian 
critical exponents may be observed [ 21 ]. 

Because of the doubts as to whether a continuum limit can be obtained, cal¬ 
culations of the spectrum of the model are sporadic. The massless photon in 
the Coulomb phase was observed in various studies [4,22-25] and the presence 
of several gauge-balls in both phases was indicated in [4]. In [23,24] the 1 
gauge-ball, the “massive photon” in the conhnement phase was found and 
in [23] its scaling investigated and found to be consistent with v = 0.33. The 
monopole mass was measured in the Coulomb phase in a dual formulation 
of the theory [26]. To our knowledge its scaling behaviour has not yet been 
studied. 


2.3 Overview of our measurements 


Motivated by the recent progress in understanding the continuum limit, in 
the present work we investigate the spectrum and its scaling behaviour. Our 
choice of 7 = — 0.2 is a compromise between the requirements of minimizing 
both the two-state signal and the autocorrelation time. 

In the simulations we have used a vectorized three-hit Metropolis algorithm 
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/3 

s^ie 

10^20 

12324 

14328 

16332 

18336 

20340 

1.100 

- 

- 

6.0 

- 

5.2 

- 

- 

1.130 

8.0 

- 

6.0 

- 

5.7 

- 

- 

1.135 

- 

- 

- 

- 

4.5 

- 

- 

1.140 

5.0 

- 

6.0 

- 

6.8 

- 

- 

1.145 

- 

- 

- 

- 

6.6 

- 

- 

1.150 

16.0 

- 

9.0 

- 

7.5 

- 

- 

1.152 

5.0 

- 

5.0 

- 

3.6 

- 

- 

1.154 

- 

- 

9.0 

- 

6.4 

- 

2.4 

1.156 

- 

- 

2.0 

- 

6.4 

- 

1.6 

1.158 

- 

- 

- 

- 

6.1 

- 

2.0 

1.159 

- 

- 

- 

- 

1.0 

- 

13.7 

1.160 

- 

- 

- 

- 

- 

- 

0.6 

1.161 

10.0 

4.0 

8.0 

1.0 

6.0 

- 

1.6 

1.162 

5.0 

4.0 

4.0 

1.0 

3.0 

1.6 

1.2 

1.165 

52.0 

4.0 

10.0 

1.0 

14.0 

0.8 

2.4 

1.170 

44.0 

22.0 

14.0 

9.4 

16.7 

2.5 

3.0 

1.180 

14.0 

4.0 

4.0 

1.0 

10.0 

0.8 

1.6 

1.190 

- 

- 

- 

- 

2.0 

- 

- 

1.200 

5.0 

- 

4.0 

- 

5.0 

- 

1.0 


Table 1 

List of our data in the confinement (upper part) and Coulomb phases. The numbers 
in the table are the numbers of measurements in multiples of 1000. 

with acceptance around 50% at every hit. The measurements have been per¬ 
formed every 25 sweeps. The used lattice sizes L^Lt, Lt = 2Ls, are listed in 
table 1. In the conhnement phase the presented data have been obtained on 
the L^Lt = 16^32 and 20^40 lattices (smaller lattices being used only for ex¬ 
plorative calculations). In the Coulomb phase the lattice size has been varied 
in a broad range. 

In table 1 we list the (3 values at which the measurements were made. The 
point (3 = 1.100 turned out to be too deep in the conhnement phase to allow 
reliable measurements of the masses. Ai j3 = 1.160 an indication of phase hips 
was observed even on our largest lattice. These two points have therefore been 
excluded from the data analysis, leaving 10 points (3 = 1.130 — 1.159 in the 
conhnement phase. In the Coulomb phase the data have been collected at 7 
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points. All data in the vicinity of (3c have been checked for absence of any 
indication of phase flips in the time evolution of the plaquette energy. 

The calculation of the GB masses at many (3 points, required for the study 
of their scaling behaviour, limited the statistics at any given (3. This made it 
difficult to distinguish GB states from background in channels with weak GB 
signal. To improve this situation we have accumulated substantially higher 
statistics for at least one point, (3 = 1.159 on the 20^40 lattice. This is the 
point closest to (3c that still has no phase flips. This allows us to determine 
the GB masses in channels with weaker signal with higher reliability at least 
at one (3. 

As listed in table 1, we have typically accumulated several thousand statis¬ 
tically independent gauge field conhgurations. We have made two sorts of 
measurements. The expectation values of the rectangular Wilson loop opera¬ 
tors Wr^t with temporal extension T and extension R in any of the spatial 
directions have been determined. All Wilson loops of sizes R = 1,..., Ls/2 and 
T = l,...,Lt/2 have been considered. From this we have obtained the static 
potential essentially by standard methods. 

Our main task, the measurement of the GB masses in various channels at 
various (3, has been performed by adopting the latest methods of glue-ball 
measurements in pure gauge lattice QGD to the U(l) gauge group. We have 
followed [27,28] in implementing the techniques of smearing and diagonaliza- 
tion of the correlation matrices. A more detailed description is given in sec. 4. 
The effective energies ej{t) in various GB channels j with zero and, in some 
cases, smallest nonzero lattice momentum have been obtained at various dis¬ 
tances t. The GB masses mj have been obtained from the plateaus of these 
energies where possible. In some channels with weak GB signal ej{t) at only 
one t, t = 1 could be used, giving at least an estimate of the GB mass. In this 
procedure the standard lattice dispersion relation has been used. 

The lattice calculations provide static potential V{R), string tension a, and 
masses rrij in the lattice units. To give these observables in physical units, e.g. 

jg necessary to specify the value and the vanishing of the 
lattice constant a. As several scenarios for the continuum limit a —0 must 
be considered, we postpone this issue to sec. 6 and until then use the lattice 
units only. 

Simultaneously with the above measurements we have also determined various 
fermionic observables in the quenched approximation. This is aimed at a study 
of the U(l) theory with staggered fermions and will be published separately 



3 Static potential 


3.1 Measurement of the statie potential 


In the static potential analysis we include Wilson loop expectation values 
whose noise to signal ratio is smaller than 20% of the expectation value itself. 
Typically we observe that Wilson loop expectation values of magnitude 0.0001 
and larger are measurable within this threshold, entirely independent of j3. 

The potential V{R) between static charges is dehned as 

r(fi)= Mm (6) 

1 —^OO J 


For hnite systems it must be specihed in a suitable way. 

In the conhnement phase the Wilson loop expectation values decay rapidly 
with increasing T. To extend the usable T interval we use the Ansatz 

{Wr,t) = , V{R) < V,{R), (7) 


at hxed R. This allows for two leading eigenvalues in the transfer matrix rep¬ 
resentation of the gauge invariant correlation functions. When compared to 
the data this truncation appears to be a reasonable approximation. The four- 
parameter fit at each R is performed solely for those values of R where at 
least 6 data points are below the 20 percent threshold in a T-interval starting 
with T = 1. We find that Vi{R) exceeds the static potential V{R) by a factor 
of about 2 throughout the conhnement phase. At /? = 1.130, the data point 
deepest in the conhnement phase, we are only able with this method to deter¬ 
mine the static potential at distances R = 1,2, 3,4. The situation gradually 
improves as the critical point is approached. At /? = 1.159, the data point clos¬ 
est to criticality, we were able to determine the static potential at distances 
1 < i? < 7. We have checked that the values of V{R) obtained in this way 
are consistent with those obtained by using the T > R loops only and setting 
Ci{R) = 0. 

In the Coulomb phase the Wilson loops of the considered sizes are much better 
measurable. There a ht to the form (7) with Ci{R) = 0 for values of T > 8 
unambiguously determines the static potential for 1 < i? < 8 on 16^32 lattice. 
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3.2 Confinement phase 


In the confinement phase we are primarily interested in the string tension a 
and its scaling behaviour. The string tension is defined through the asymptotic 
behaviour of the potential V{R) oc aR as R ^ oo. In order to extract a at 
the available distances R, a reliable parameterization of the short distance 
part of the potential would be of much help, similar to lattice QCD. In the 
U(l) theory we face a serious obstacle: the lack of theoretical knowledge of 
the potential at short distances in the strongly interacting theory which is not 
asymptotically free. 

We have experimented with various fits to the potential on the whole available 
/^-interval, parameterizing the short distance behaviour of the potential in 
several ways. The lattice Coulomb potential of a massless exchange particle, 
as well as the lattice Yukawa potential of a massive exchange particle were 
considered. At the available distances these Ansatze influence the values of a 
significantly. This effect is at the 10 percent level for our smallest /5-values, 
where the potential is only known at few i?-values, and is less pronounced, 
at the level of a few percent, close to the critical point. This is so also if 
the gradient of the potential, the force, is considered. In summary, we cannot 
avoid a dependence of the fitted string tension values on the Ansatz for the 
short range part of the potential if all the data for V{R) are fitted. As this 
systematic error decreases when the critical point is approached, it distorts 
the scaling behaviour of a determined in such a way. 

In this situation we expect that a straight line fit to the potential for the largest 
i?-values results in a less ambiguous and model independent determination 
of the string tension. For the final determination of a at all the values 
considered we used the potential data at the 3 largest R values with V {R) 
being measurable in the sense described above. In table 2 the obtained values 
of the string tension in the confinement phase are listed. 

The square root of the string tension y/a is expected to scale like a mass. 
We fit the string tension data using the critical value fic = 1-1607 which is 
obtained below from the scaling of the GB spectrum. Recalling that our string 
tension determination is less affected by systematic uncertainties close to the 
phase transition we restrict the fit to 1.14 < /5 < fic- The fit to the scaling 
form 




( 8 ) 


results in the exponent value = 0.35(2) and a x^/Ydf- value for the fit of 
about 1.9. We test the stability of the fit by the omission of further data points 
and obtain //-values scattering in a 10 percent interval around the quoted value. 
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p 

a 

16332 20^40 

1.130 

0.268(6) 

1.135 

0.244(5) 

1.140 

0.226(2) 

1.145 

0.178(4) 

1.150 

0.137(6) 

1.152 

0.109(10) 

1.154 

0.118(11) 0.139(15) 

1.156 

0.086(14) 0.082(13) 

1.158 

0.066(9) 0.062(17) 

1.159 

0.028(13) 


Table 2 

Values of the string tension a in the confinement phase. 

The square root of the string tension data and the fit are displayed in fig. la. 
Within the present numerical and systematic uncertainties we conclude that 
the string tensions scaling behaviour is consistently described by the non- 
Gaussian exponent = v^g- 

3.3 Coulomb phase 


The static potential in the Coulomb phase is well known to be dominated by 
the exchange of the massless photon. Its determination and a confirmation 
of its Coulomb form is a spin-off of our study of the spectrum in this phase. 
Nevertheless, it is of interest to determine the maximal renormalized coupling 
at 7 = —0.2, as this has not been done before. 


An appropriate representation of the lattice Coulomb potential Vl,c is ob¬ 
tained by the lattice Fourier transform of a massless bosonic propagator. It 
has the form 


Vl,c{R) 


Att _ plkiR cy 


,9) 

Ls 


We have htted the static potential at seven /^-points in the Coulomb phase 
with the form 


V{R) = -aKVi^,c{R) + (TR + C ( 10 ) 
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p 

LpLt an 

1.161 

20^40 0.180(5) 

1.162 

20^40 0.186(3) 

1.165 

15332 0.179(4) 

1.170 

15332 0.169(3) 

1.180 

15332 0.162(2) 

1.190 

15332 0.157(2) 

1.200 

15332 0.152(2) 


Table 3 

Values of the renormalized fine structure constant an in the Coulomb phase. 

and obtain a good description of the potential. The coefficient a is consistent 
with zero within error bars, conhrming the vanishing of the string tension 
in the Coulomb phase. We therefore omit the linear term altogether in the 
determination of the renormalized hne structure constant q;r. The values of 
ctR grow when the critical point is approached. Our data are collected in table 
3. 

Based on the analogy between the low temperature (small coupling) expan¬ 
sions of the D = 2 XY spin model in its spin-wave phase, and of the compact 
U(l) gauge theory in the Coulomb phase, one expects [14] that the renormal¬ 
ized fine structure constant behaves according to 

an{/3) = q ; r,c - Ac,{/3 - (3^)^ , p > pc- ( 11 ) 


At the critical point a hnite value of ctr^c has been predicted and it has been 
suggested that this value might be universal [13,14]. This is in complete anal¬ 
ogy to the D = 2 XY model, where the corresponding coupling of the 2- 
dimensional analog of the Coulomb potential has the critical value 1/4. 

Previous simulations [15] of the model with the Wilson and Villain actions in 
the Coulomb phase showed consistency between the data and (11): for both 
actions obtained the value ctr^c — 0.20 with an error about 0.02. 

In order to determine q;r,c we extrapolate the values of q;r obtained at non- 
critical P values in the Coulomb phase to Pc via (11). We have used the value 
for the critical coupling Pc, as determined by the scaling of the gauge ball 
spectrum. With three free parameters aR,c, Aa and A we have obtained 

aR,c = 0.19(1) , A = 0.5(2) (12) 

with x^/-^DF = 0.9. We hnd that q;r,c is quite well determined by this extrap- 
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olation within the Coulomb phase (cf. £g. lb). 

An alternative determination of the renormalized fine structure constant at 
Pc can be obtained by assuming that the short range properties of the static 
potential are smooth at the critical point. Thus in the vicinity of the critical 
point it is sensible to dehne q;r from the short range part of the static potential 
in the confinement phase too; q;r,c can then be obtained by interpolation. This 
gives results in agreement with ( 12 ). 

Our result for ctr^c at 7 = —0.2 is consistent with those obtained for other 
actions [15] and thus further supports the conjectured universality of q;r,c in 
the pure U(l) gauge theory. 




Fig. 1. (a) a in the confinement phase. The line is the fit ( 8 ) of data at /3 > 1.140. 
(b) Or in the Coulomb phase. The line is the fit (11). The lattice sizes are 16^32 
(squares) and 20^40 (crosses). 


4 The gauge-ball measurements 

4.1 Quantum numbers 


Our GB operators for a given time-slice are composed of appropriately sym¬ 
metrized loops containing links only in the three spatial directions. The com¬ 
binations of loops must have a well-defined behaviour under the point group 
for cubic symmetry with inversion, Oh- For rotations one has the proper cu¬ 
bic group O with its five irreducible representations R = Ai, A 2 , E, Ti and 
T 2 of dimensions 1, 1, 2, 3 and 3 respectively. Inversion then introduces the 
parity quantum number P = ± 1 , involving direct group product such that Oh 
behaves as O ® P. 

As the representations are complex, we can take the real or imaginary part 
of the loop, corresponding to charge parity C = -|-1 and —1, respectively. 
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irr. rep. 

dimension 

smallest spins 


1 

0, 4, ... 

^2 

1 

3, 6, ... 

E 

2 

2, 4, ... 

Ti 

3 

1, 3, ... 

T2 

3 

2, 3, ... 


Table 4 

The smallest spins that are contained in the irreducible representations of the cubic 
group [30]. 


type 

^ links 

^ orientations 

shape 

plaquette 

4 

3 

EJ 

rectangle 

6 

6 

/ / / 

chair 

6 

12 


flange 

8 

48 

^. 7 


Table 5 

Loops used for the construction of the gauge-ball operators. Their shapes are ex¬ 
tended only in the spatial directions. 

Thus a complete categorization of a state is with the hve possible R as 
above [30]. 

In the following we consider all these 20 states with momentum p = 0 as well 
as some with p = 1 (in units of 2t[/L g). We denote these channels by 

3 = (13) 


The representations R contribute in the continuum limit to spin as indicated 
in table 4. In glue-ball calculations one tends to identify the lowest energy 
lattice state with the lowest continuum spin to which it contributes, but this 
is not necessarily correct. 

For p = 1, the states no longer have exactly the symmetry of Oh represen¬ 
tations as there is no longer a dehnite parity. Thus in general our non-zero 
momentum states have been chosen to minimize the problem of mixing be¬ 
tween states. In most channels the p = 1 results simply conhrm those for 
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p = 0. In these cases either the mixing is negligible or the mixed states have 
higher masses than the ground state of the ‘native’ channel. Two exceptions in 
the confinement phase and one (the photon) in the Coulomb phase are noted 
in our discussion of the results. 

We have used the loops shown in table 5. In particular the eight-link operator 
(flange) has no remaining symmetry and so contributes to all . One can 
calculate how the various loop orientations project onto states of Oh using 
the tables in [30,31]. Thus for zero momentum states we have used suitable 
combinations of loops to represent the representations of Oh and determined 
the lowest lying states in all channels. In order to identify the massless photon 
in the Coulomb phase we have to use a non-zero momentum. We discuss this 
further below. 


4-2 Improved operators and effective energies 


To improve the operators we have used the standard factor-of-two blocking 
prescription [27]. Accordingly, sets of links on the original lattice are combined 
together to form links of twice the length by adding the contributions from 
the central link of length 2a plus all four spatial “staples” of length 4a: 


{x, pI) = U{x, p)U{x + fi, n) 


y~] U{x, i')U{x + n, p)U{x+O+p, p)U\x + 2p, u) 


Spatial 

+W{x — u, v)U{x — t)^ ia)U{x — u+fi, jji)U{x — t)+2fi^ n) 


(14) 


The resulting complex number is normalized back to an element of U(l) with 
unit modulus. This operation is essentially a block-spin renormalization trans¬ 
formation to a coarser lattice and can be performed recursively, generating 
f/(2) from and so on. The practical limit is reached when the Wilson loop 
combinations constructed from these blocked links reach the size of the spatial 
lattice. 


This gives us many operators 0”(t) for a GB state j at Euclidean time t, 
where the index r now labels all contributions from different blocking levels 
and loop shapes to a given j. From this one can form a connected hermitean 
correlation matrix, 

C7(t) =i:i{or(«o+«)o*(«„)> - {o"((„+«))(oj((„))]. (15) 

^0 


The correlation matrix (at least for a complete set of states) allows for a 
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spectral decomposition 


q^(t) = ^ (e-®^ + , (16) 

k 

where we assume nondegenerate eigenvalues. The coefficients represent 
the projections of our GB operators on the energy eigenstates. Since we have 
only a limited number of states we have to assume that the diagonalization of 
the truncated matrix still gives the physical energies to a good approximation. 
The dominant contribution to the matrix at large t comes from low mass states, 
which are the ones we are most interested in. 

We now diagonalize the correlation matrix by analogy with [28,32,33]. The 
actual operation we perform can be written as a generalized eigenvalue prob¬ 
lem. 


Cj(to + l)uj — \jCj(to)vj 


(17) 


in which we pick the eigenvector Vj corresponding to the smallest energy, 
and hence the largest eigenvalue \j (the rs indices have been suppressed). 
Solving this is simplihed by the fact that the Cj{t) have all been symmetrized. 
The larger we choose to here, the smaller the contribution from higher states. 
However, the statistical fluctuations increase rapidly with to and in practice 
the gain in stability by picking to = 0 is more important when looking only at 
the lowest few states. The corresponding Vj determines the linear combination 
of operators used for all subsequent analysis at every t. Its correlator Cj(t) 
gives the effective energy and the full correlation matrix is not required again. 

Another practical point is that many of the contributions to Cj actually have 
very similar sets of coefficients A^^ for these lowest states, making the diago¬ 
nalization less stable. It should be remembered that the index j labels blocking 
level as well as the Wilson loop shape: we have reduced the number of opera¬ 
tors which appear in the correlation matrix by picking out the loop shape with 
the smallest 6^(0) = log(GJ''(0)/GJ''(l)) from each blocking level, and simply 
discarding the rest. Thus between 3 and 5 states actually remain in the cor¬ 
relation matrix, depending on the lattice size. The effect of the blocking — 
the projection onto the lowest state increases (ej(0) decreases) and then falls 
off again as the amount of blocking is increased — means that this provides 
a good spread of Ay^ and the diagonalization in all cases works well. Where 
a comparison is possible the results for the lowest mass are indistinguishable 
from those with the full correlation matrix. We have found that the blocking 
and diagonalization of the correlation matrix works essentially as well as has 
been found in pure gauge QCD [28]. 

The effective energies ej{t) are obtained from the correlation Cj{t) by solving 
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numerically the following two coupled equations for ejit) and Kj{t), 

c,(t) = (18) 

Cj{t + 1 ) = 

Their approximate solution for small t is 

4.3 Errors of effective energies 


The errors of the effective energies ej{f) were obtained directly by the boot¬ 
strap method. In this procedure, data consisting of sets is sampled to 
produce a new ensemble of Nd sets; each sample is completely random, so 
that the new ensemble can contain any of the original sets more than once 
or not at all. This ensemble is then processed in exactly the same way as the 
original. The procedure is then repeated an arbitrary but suitably large num¬ 
ber of times iV^^, producing a spread of results. The standard deviation 
of these results gives a reliable estimate of the statistical error in the original 
ensemble. 

Our initial data was blocked down by a factor 50, producing an Nd of typically 
a few tens of data sets for each observable. We have picked Ws = 99- We have 
checked that using the same variational basis for all the bootstrap data, rather 
than re-diagonalizing the correlation matrix for each bootstrap sample, does 
not affect the errors obtained, and have therefore performed the bootstrap on 
the data only after the diagonalization. 

We have arbitrarily chosen only to look at the hrst eight correlations and 
hence seven ej{t) for t = 0 to 6. The choice of the range of t over which we 
perform hts to extract the real energies Ej is described in the next section. 

We have analyzed the appropriate range of effective energies in two ways: 
hrstly, using a one-parameter uncorrelated £t to the data over the chosen 
range of f, and secondly by taking a weighted mean of all the values over 
the range. Using the second method we processed the bootstrap samples with 
the same weights, giving an estimate of the overall statistical error which is 
presumably more reliable than that on the fit parameter. Certainly the errors 
with this method appear consistently larger than with the other; some effective 
energies have an error larger by a factor of around two. This is particularly 
noticeable in the states near to the phase transition with a good signal. We 
have therefore used the values from the bootstrap of the weighted means for 
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our statistical errors. 


This analysis gives essentially statistical errors only. An estimate of the sys¬ 
tematic errors of our hnal results will be made at the end of the next section. 


5 Scaling of the GB masses in the confinement phase 


5.1 Gauge-ball energies and masses 


The list of GB channels we have investigated is given in the first column of 
table 6. The momenta considered are indicated in parentheses. In the second 
column we give the quantum numbers of the expected continuum state with 
the smallest spin, in accordance with the table 4. 

In two cases our results for the conhned phase should be interpreted as showing 
a clear signal for mixing when the momentum is non-zero, and hence the naive 
interpretation is not correct. These states, the and T 2 , are marked 
with an asterisk in table 6. For p = 0 they give no clear signal, increasing the 
likelihood that states mixed from other will be visible. To see that the 
mixing is as expected, one needs to consider the behaviour of Oh states when 
the symmetry is broken along a particular axis, as by our (1,0,0) x 27r/Ls 
momentum boost. The reduced symmetry group is dihedral: one finds that 
the component whose own axis of symmetry lies along the momentum 
axis (i.e. longitudinally polarized) has the same behaviour under the reduced 
symmetry group as the Af^, and that the T 2 likewise mixes with the T 2 "*'“. 
This agrees with our results. Also the channel T 2 "’''(l) might have an admixture 
from one of the 2 ^^ states (see below). 

In each channel we have determined the effective energies ej{t) for as many 
t < 6 as possible. In principle one should obtain each GB energy Ej by a £t 
to their plateau which is expected at large t. In practice, in some channels 
plateaus can be found only at moderate f, since ej{t) has very large errors or 
is unmeasurable for larger t. The quality of the plateaus is therefore different 
in different channels and varies with (3, being best close to fdc- Sometimes we 
cannot identify a plateau at all. 

We therefore classify the channels according to the quality of the signal for a 
definite GB energy Ej, and extract this energy in somewhat different ways. 
The classification is indicated in the third column of table 6; this holds only 
for the confinement phase. The meaning is as follows: 

states: At most f3 a plateau in ej{t) was found and fitted at t > 2. 
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Table 6 

GB states with momentum p = 0,1 (in the units LsI2tt) observed in various channels 
in the confinement phase and their scaling behaviour. If two momenta are listed 
in one line, the mass has been assumed to be the same. The quantum numbers 
of the continuum states with smallest possible spin are indicated in the column 
“continuum”. The clarity of the effective energy plateau and the reliability of the 
mass determination is indicated by a corresponding number of “+”. In the channels 
denoted “?” the evidence for a mass is not clear. In the columns ‘T” and “mass” 
the scaling exponent of the GB mass and its approximate proportionality to one of 
the mass scales is indicated. The ratio of each mass to the corresponding mass scale 
is given in the column “ratio”. For the explanation of the continuum states marked 
with * see the text._ 



continuum 

quality 

u 

mass 

ratio 

yl++(0,l) 

0++ 




1 

Tf+(1) * 

0++ * 



mg 

1.01(6) 

^r(o) 

3+- 




1.00(4) 

E+-i0) 

2+- 

+ + + 



0.98(4) 

r+-(o,i) 

1+- 



mng 

1 

TtiOA) 

2+- 




1.00(4) 

T,-{1) * 

2+- * 




1.04(5) 

E++i0,l) 

2++ 




2.06(9) 

r2++(0,i) 

2++ 

++ 


~ 2 mng 

2.07(9) 

ylr+(0) 

0-+ 




2.88(13) 

yir-(O) 

0— 


Z/ng 


2.73(14) 


3++ 




3.22(16) 

^2”“(0) 

3— 




2.73(15) 

E-{0) 

2— 




2.95(14) 

r++(o,i) 

1++ 

+ 


(2.6 — 3.6)mng 

3.61(18) 

Tf-(0,1) 

1— 




2.94(13) 

r2-+(o) 

2-+ 




3.10(14) 

r2-+(i) * 

2++ * 




2.63(12) 

r2-(o) 

2— 




3.25(16) 





1 / 3mng 

0.34(2) 

yl+-(0) 

0+- 




^2^(0) 

3“+ 

? 

? 


E-+{0) 

2“+ 




rf+(o) 

1-+ 
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++ states: At most (3 a plateau in ej{t) was found only if t = 1 was included. 
This plateau was then fitted. 

+ states: A plateau in ej{t) was found only at the highest statistics point 
f3 = 1.159 if t = 1 was included. This plateau was then htted. For other f3 
the value of ej{l) was taken for the GB energy Ej. 

? states: No plateau in ej{t) was found at t > 1 even for (3 = 1.159. Then we 
do not obtain any GB energy. 

In order to illustrate the quality of the data, we show in fig. 2a results for the 
effective energies of (+ + + quality) and (++ quality) with p = 0 
obtained at a typical data point (3 = 1.145 on the 16^32 lattice. The effective 
energies of the same states at the data point with best statistics and smallest 
errors, (3 = 1.159 (close to j3c) on the 20^40 lattice, are shown in £g. 2b. The 
GB energies Ej with errors are shown as horizontal lines. 



0 1 2 3 4 5 6 



0 1 2 3 4 5 6 


t t 

Fig. 2. Effective energies from A^^(O) and T2~^{0) correlation functions (a) at 
P = 1.145 on the 16^32 lattice and (b) at /3 = 1.159 on the 20^40 lattice. The 
plateaus were fitted omitting the points at t = 0,1 for and the point t = 0 for 
The GB energies Ej and their errors are indicated by the horizontal lines. 


The GB energies Ej at p = 0 and p = 1 have been used to calculate the GB 
masses rrij by the lattice dispersion relation 


3 

2{cosh. Ej — 1) = + 2^(1 — cosp^). (20) 

ti=i 


If the masses found in the p = 0 and p = 1 channels with the same 
turned out to be consistent, we have continued the analysis assuming their 
equality, i.e. they were fitted by exactly the same parameters. In the table 6 
these “pair” channels are represented by one line only, as well as by the same 
symbol in the hgnres 3 and 6. 
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Fig. 3. Values of Vj obtained in each channel separately. The double vertical line 
separates two groups with distinctly different Vj. The dashed vertical lines separate 
channels with different quality signal within the group with non-Gaussian exponent 
values. The ordering within this group is with decreasing reliability (table 6) from 
left to right. 

5.2 Two mass scales 


Having determined the GB masses for various (3 we have investigated their 
scaling behaviour with j3. First the fits of the form 

= ( 21 ) 


were performed for each GB channel j individually. We found two groups of 
masses with strikingly different scaling behaviour. Within the whole (3 range a 
large group of the GB masses scales with roughly the same exponents Vj close 
to the non-Gaussian value 0.365(8) found in [8-10]. However, in a small group 
consisting only of the , p = 0,1 and the , p = 1 channels the values 
of Vj are approximately Gaussian, i.e. 1/2. The values /?/ are quite consistent 
with each other in all channels. Therefore, in the further analysis we have 
assumed the same value of (3c in all channels. 
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0.726 +- 0.022 


m(A/") 


Fig. 4. A log-log plot of the masses of the and gauge-balls in the interval 
1.13 < f3 < 1.159. The results from the channels with p = 0,1 are denoted by 
circles and squares, respectively. The solid line is the result of a straight line fit 
considering errors in both directions. The dashed straight line with a gradient of 
one, i.e. assuming equal exponents, is shown for comparison. 

The results for Uj obtained from (21) with common {3c are shown in fig. 3. A 
clustering around two values is clearly seen. In order to further demonstrate 
the difference of the scaling behaviour between the two groups we plot the 
masses of the clearest members of each group, the T^~ and channels 
with p = 0,1, in a log-log plot against each other in £g. 4. The data have a 
slope distinctly different from one. 

Naturally, as seen in fig. 3, the errors of uj in channels with less accurately 
determined masses are rather large. However, for example in the T^~ and 
channels with p = 0 and p = 1 the masses are quite accurate and the 
statistical errors are small. Performing fits to each of these two channels with 
both momenta we obtain 


^(T+“,p=o,i) ~ 0.37(3), 

= ( 22 ) 

The masses in these two channels with both momenta (distinguished by dif¬ 
ferent symbols) are shown in fig. 5. The corresponding fits by means of (21) 
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Fig. 5. Masses of the (lower points) and glueballs in the confinement 

phase versus (3 for both momenta (p = 0,1 being denoted by circles and squares, 
respectively). The dashed curves are fits to the data shown. The full curves represent 
the mass scales (26). 

are indicated by dashed lines. 

These observations suggested the next step of the analysis: we have assumed 
the same exponent Vj for each group, denoted and z/g for the non-Gaussian 
and Gaussian group, respectively. Then a joint £t was performed with these 
two exponents, common Pc, and the individual amplitudes Cj as free parame¬ 
ters: 


mj = CjT''', f=iig, g. 


(23) 


The resulting values of the exponents are 
z/ng = 0.367(14), 

z/g = 0.51(3). (24) 

They are essentially determined by the most accurately determined masses, 
such as those in the T]^“(0,1) and Aj'“’''(0,1) channels. The critical point at 
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7 = —0.2 determined in this way is 


/3c = 1.1609(2). 


(25) 


The amplitudes, giving the ratios of the masses in each group, are described 
below. 

All the reliably determined GB masses are well described by this £t. Those 
which are not fully reliable are at least consistent with it. From this we con¬ 
clude that the system has two mass scales which we denote by rUng and mg. 
Each GB mass scales according to one of these mass scales. We choose the 
fitted scaling behaviour of the GB operators Ti~ and to define the two 
mass scales: 


^ng — CngT Cng — C(T+“,p=0,l) ~ 5.4(3), 

mg = Cg = = 7.4(6). (26) 


The full curves in fig. 5 represent these mass scales mng and mg. We observe 
that the dashed curves are almost indistinguishable from the mass scales (26). 
Our determination of these scales is thus not sensitive to the selection of only 
one channel (with both p) data for their definition, nor to the use of a joint 
value for Pc- 


5.3 Gauge-ball spectrum 


Table 6 summarizes the results for the GB masses in the confinement phase 
and gives the observed masses in multiples of the two “standard” masses mng 
and mg. For the states denoted by -|- our assignment is only tentative. 

The results for the individual masses are presented in fig. 6 and in the last 
two columns of table 6. In fig. 6 the individual amplitudes Cj are shown. 
They summarize our results for the masses obtained from all P points in 
the confinement phase and reflect their scaling behaviour. It is apparent that 
the states in the non-Gaussian group cluster around three values of Cj, being 
multiples of the lowest mass states in this group. 

An alternative way to determine the GB masses is to rely only on the high 
statistics data point, P = 1.159 on the 20^40 lattice. The absolute values of 
rrij at this point are presented in fig. 7. Here we show the results for p = 0 and 
p = 1 separately in all channels, including the pairs. Of course, at one fixed 
P the difference in the scaling behaviour is not observable. However, keeping 
in mind that there are two groups of states, one can compare figs. 6 and 7 for 
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Fig. 6. The amplitudes Cj of the fits (23). The vertical lines have the same meaning 
as in fig. 3. 

each group separately and find an agreement between the mass ratios. The 
two-fold and perhaps three-fold clustering of the values in the non-Gaussian 
group around the lowest mass is again apparent. 

As seen in fig. 7, the dispersion relation (20) is slightly violated for the pairs 
i7^’''(0,1) and T2^’'"(0,1). Motivated by the fact that the corresponding masses 
are very close to 2mng, we have also tried to apply the dispersion relation 
assuming that it is really two-particle states which are seen in these channels, 
either both particles having zero momentum or one of the particles moving 
with the lowest nonzero momentum. Indeed, the £t of Ej in these channels 
gave the mass rUng for each of the two particles even for p = 1. This result 
is indicated by crosses in fig. 7. It suggests that the i7’''+(0,1) and T2*'^(0,1) 
channels are dominated by states of two gauge-balls from the non-Gaussian 
group with mass rring. Of course, further verification of such a possibility is 
necessary. 

Only the ratios of the amplitudes within each group are relevant for the con¬ 
tinuum limit, and therefore we show in the last column of table 6 the ratios 
Tj to the corresponding mass scales (26), 

nij = TjUii, f = ng, g. (27) 
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Fig. 7. GB masses rrij at the point with the highest statistics, /3 = 1.159 on the 
20^40 lattice. The p = 0,1 results are indicated by the circles and squares, respec¬ 
tively. The crosses belonging to the 1) and r^^(0,1) channels are masses of 

individual particles when the energies Ej in these channels with p = 1 are assumed 
to correspond to two-particle states with lowest nonzero momentum. The vertical 
lines have the same meaning as in hg. 3. 

Relating the different scales set by the string tension and by the mass rUng 

= Tuning, (28) 

we obtain 

r, = 0.34(2). (29) 

This value is given in table 6 for a comparison with the GB masses. 

It is interesting to compare our results in the confined phase with those of 
quenched QCD (i.e. pure gauge SU(3)), for which the most recent analysis 
of all appeared in ref. [34] (see in particular £g. 1), and of the and 
also in [35]. Although our results here are spread over a larger variety of 
lattice sizes and couplings, the higher efficiency of U(l) calculations implies 
that the statistical accuracy of the results is comparable in the two cases. 
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As in SU(3), we find that masses tend to appear together in bands. This is 
more pronounced in our case as it appears in low-lying and therefore well 
measured states. In both cases we hnd the lowest lying state to be the most 
symmetric, namely We also hnd a low-lying 2'^^ (lattice and 
state. A clear difference to SU(3) is that in the U(l) case there is a band of 
intervening PC = H— states for spins 1, 2 and 3, all with similar mass, which 
do not appear in SU(3) (though in the latter case there is a state with 
approximately twice the lowest glueball mass). Indeed, these states in U(l) 
are the lowest ones with the non-Gaussian scaling behaviour. On the other 
hand, we do not hnd the low-lying O"’*' (and maybe 0^ , which had a large 
error) of SU(3). 

Heavier states have correspondingly large errors, making any further compar¬ 
ison too difficult. It is at least clear that the similarity of conhning behaviour 
is not by itself enough to guarantee a similar spectrum. 

One further similarity to QCD is remarkable: It has been known for long that 
in the lattice QCD the (3 dependence of the 0'''+ glueball is somewhat diherent 
from that of the other glueball states. We shall return to this point in the next 
section. 

The simplest interpretation of the observed spectrum might be the following: 
The Gaussian spectrum consist of one scalar only. The non-Gaussian 
group has a nontrivial but quite simple spectrum consisting of the , 2^ 
and d"*"" states of the same mass rUng. Higher spin assignments for the J^~ 
states are not ruled out, however. In the other channels the masses are ~ 
2mng and (2.6 — 3.6)mng. They may well correspond to two- and three-particle 
states of the lightest non-Gaussian gauge-balls. The multi-particle states of the 
Gaussian scalar do not seem to contribute signihcantly to the non-Gaussian 
channels in spite of its low mass. This further supports the expectation that 
the belongs to a “trivial” theory. 


5-4 Error analysis 


The errors given up to now in this section have been obtained by purely 
statistical procedures. They are based on the essentially statistical errors of 
the energies Ej and the corresponding masses ruj, described in the previous 
section. As seen in (24) and (25), they are rather low, on the level of 4% for 

Wg- 

However, it is obvious that various additional errors - partly of systematic 
origin - should be expected. The main reason is that we have to stay away 
from the critical point in order to avoid hnite size effects and the two-state 
signal. The value of Pc then enters very sensitively into the determination of 
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V by fits like (23). Furthermore, mass values rrij ~ 1 have to be used in such 
a fit. This raises the question of nonleading terms in the scaling behaviour. 

To get at least an estimate of these effects, we made several analyses of the 
scaling behaviour of the well determined (+++) masses when jSc is approached, 
modifying the procedure in various ways. In most cases the obtained values 
for u and Pc varied slightly with the procedure. 

(1) We modihed the P interval in which the power law hts were performed, 
dropping up to hve P points most distant from Pc- No signihcant change 
of the //-values has been observed as long as the same Pc was used. 

(2) Instead of (20) we have used the dispersion relation with replaced 
by 2(coshmj — 1) (this amounts to ruj = Ej for p = 0). The difference 
between these two frequently used relations is on the O(m^) level and 
thus indicates a possible effect of nonleading terms. 

(3) We analyzed separately the masses from the p = 0 and p = 1 channels. 

(4) We analyzed separately the masses obtained on the 16^32 and 20^40 
lattices. This also corresponds to a variation of the P interval close to Pc- 

(5) The same data were fitted both with free Pc and with this value fixed 
within the interval of its values obtained by other modihcations. 

( 6 ) The Gaussian and non-Gaussian group was htted with independent Pc 
for each group. The difference of Pc values is about or less than 0.0002. 

Dropping the extremal valnes for each parameter, we found their valnes lay 
in the following intervals: = 0.32-0.38, z/g = 0.41-0.57 and Pc = 1.1605- 

1.1609. The ratio ~ 0.62-0.80 is always clearly different from one. These 

intervals are signihcantly broader than the statistical uncertainties. Assnming 
some peaked distribution of the valnes in these intervals, we take one quarter 
of their widths as an estimate of the systematic errors not acconnted for by 
the purely statistical analysis. 

The valnes of the critical exponent v given in (1), (2), (4), and of Pc given in 
(3) are the central points of these intervals. The qnoted errors are the simple 
sums of the above systematic errors and typical statistical errors. 


6 Three scenarios 


The occnrrence of two different correlation length exponents may seem surpris¬ 
ing in the light of the simulations on spherical lattices [ 8 - 10 ]; there the scaling 
of several bulk observables could be well described by means of only one expo¬ 
nent, z/ng. In particular, the Fisher zeros at 7 = —0.2 showed no deviation from 
the asymptotic hnite size scaling behaviour determined by this value of the 
exponent z/. However, that analysis was based completely on the bnlk energy. 
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In this observable there is little chance of hnding nonleading behavionr. The 
contribntions scaling with smaller critical exponents will be dominant near 
the phase transition. In terms of the eigenvalnes Aj of the linearized renormal¬ 
ization gronp matrix this is qnite nnderstandable: since A* = (where s is 
the scale change factor), the smallest vi dominates. 

Therefore we see no contradiction with the earlier resnlts in the appearance 
of two different correlation length exponents in onr present simulations. Fur¬ 
thermore, it does not appear unnatural when the possibility of the presence 
of a TCP at 7 = 7 o is taken into account. The TCP’s are known to have 
exponents in general different from those associated with the adjacent ordi¬ 
nary critical lines [17], and they tend to dominate the scaling behaviour in 
their vicinity when the critical manifolds are approached from outside. The 
cross-over regions between tricritical and critical scaling behaviour may ex¬ 
tend quite far away at a small angular distance from the critical lines [36,17]. 
Furthermore, different observables, e.g. correlation lengths, may have different 
crossover regions.^ 

However, when trying to explain the observed scaling of the GB masses at 7 = 
—0.2 with help of the TCP, we see at present no reliable possibility of saying 
which one, if any, of the two observed n-values corresponds to the tricritical 
scaling. Instead of guessing we formulate three scenarios, which hopefully can 
be tested in future simulations. In each of them we assume that 70 > —0.2. 

In all the following scenarios one can consider two different continuum limits, 
depending on the mass = mj/a that one chooses to £x in physical units. 

If it is one of the non-Gaussian group, then the lattice constant vanishes as 

«(^) = ^ ( 30 ) 


and analogously for the Gaussian group. 


a(r) 


nip 


phys 

mg 


oc 


(31) 


Scenario C: 

The point 7 = —0.2 is a critical point with two relevant mass scales, i.e. 
both exponents and z/g are critical. Two different scaling laws at the same 
critical point are unusual, but not, to our knowledge, impossible. A somewhat 
analogous situation is known in the three-dimensional pure compact U(l) 


^ We thank D.P. Landau for a helpful advice in this question. 
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gauge theory [37]. As the ratio of masses approaches zero, 


TTlng 


oc r 


oc r 


0 . 5 - 0.365 


0 , 


(32) 


we can consider two different continuum limits: 

(1) Cng: Keeping constant, the spectrum consists of those GB states 

which scale with z/ng as states with hnite nonzero mass, and the string 
tension is hnite. The GB would be present as a massless state. This 
would be an interesting continuum theory. 

(2) Cgi Keeping constant, only the O’*"’'' state can have a hnite mass, 

whereas the masses of all other states and the string tension would run to 
inhnity and decouple. Presumably this would be a Gaussian theory with 
a noninteracting scalar. 

This scenario does not make any use of the special properties of a TGP and 
might hold if there is no TGP with special tricritical exponents at j = jq. 

In the tricritical scenarios discussed next, we assume the tricritical point at 
7 = 7 o and take into account the possibility that some of our scaling results 
are not yet asymptotic, but obtained in a domain where a precocious form of 
scaling with a “wrong” scaling exponent is obtained. 


Scenario Tgi 

The tricritical point is Gaussian, i.e. z/g is tricritical and z/^g is critical. The 
continuum limit at 7 = — 0.2 is like that in the scenario Cng except for the mass 
of the GB state. That mass shows “wrong” (tricritical, i.e. nonasymptotic 
at this 7 ) scaling behaviour for those distances from [3c we were able to inves¬ 
tigate on lattices of limited size. In this scenario it would have to change its 
scaling behaviour closer to the phase transition, adopting the exponent z/ng. 
This state could thus have a nonvanishing mass. The data close to (3c suggest 
that this mass would be small relatively to mng. 

At the TGP the continuum limit would be Gaussian and our present results 
would say nothing about the spectrum there except the presence of the light 
0++ GB state. 


Scenario Tng: 

The tricritical point is non-Gaussian, i.e. z/ng is tricritical and z/g is critical. 
At 7 = —0.2 all the GB masses except that of the O’*"’*' state, as well as the 
string tension, show “wrong” (tricritical) scaling behaviour on our lattices but 
would change to critical scaling when approaching the phase transition more 
closely. All masses and the string tension stay hnite in the continuum limit. 
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However, because of the Gaussian value z/g the theory might be trivial with 
logarithmic corrections which we cannot detect. 

A nontrivial continuum theory analogous to Cng would then be expected at 
the TCP. As in Tng at 7 = — 0.2 the state might have a small mass. 


At the moment it is difficult to decide which of these scenarios is the correct 
one, but the most important conclusion holds in any of them: a continuum 
limit is possible in which the observables scaling like rUng remain hnite and 
nonzero in physical units. Their ratios are given by the factors r in front of 
m^g in eqs. (27) and (28), and in table 6 . This nontrivial limit is obtained on 
the critical line in the scenarios Cng and Tg, whereas in the scenario Tng it is 
obtained at the hypothetical TCP. The ambiguity which remains is the mass 
in physical units of the O’*"’'' GB state in such a continuum limit: it is zero in 
the Cng scenario whereas in both the T scenarios it may be nonzero but small 
relatively to m^g. 

Both scenarios T imply that in the range of correlation lengths we were able to 
investigate one of the two groups of channels shows a behaviour at 7 = — 0.2 
which actually corresponds to the TCP, and not to the critical point at this 
7 . This can be explained assuming a broad dominance angle of the TCP for 
the corresponding states. However, we hnd it remarkable (peculiar?) that the 
would-be critical value of (3 of the group with “false” tricritical behaviour 
agrees within our precision with that of the other group. 

The special behaviour of the (strictly speaking A^’*') GB suggests a possi¬ 
ble parallel with pure SU(iV) gauge theories with a mixed action. In that case 
the Wilson action, containing the trace of the plaquette in the fundamental 
representation, is supplemented by a term where the trace of the plaquette is 
taken in the adjoint representation with coupling 7 . The procedure is analo¬ 
gous to the U(l) with extended Wilson action. For SU(2) and SU(3) the theory 
has a hrst order line with an endpoint at positive 7 . In the SU(3) case it was 
found [38] in this region that, while the mass remains hnite, the 
mass goes to zero. The critical exponent is not known. Other masses have not 
been calculated, although it is a long-standing result for SU(2) [39] that the 
string tension also remains hnite around the corresponding end point. Thus 
the seems to be singled out for special behaviour as it is in the results 
presented here. In SU(3) there is presumably no phase transition for values of 
7 below the endpoint, apart from a near-coincidence with hnite temperature 
ehects on small lattices which appears to be accidental [40][3]. It is therefore 
natural to suppose the special behaviour is due to the endpoint itself. There 

^ The possibility of a different interpretation, very similar to the phase diagram of 
the U(l) theory, has been pointed out in [41]. 
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are other examples in lattice models where such endpoints are associated with 
Gaussian behaviour (see e.g. [42]). Depending on the choice from the above 
scenarios where the gauge-ball in the U(l) theory scales, either the TCP 
or the critical line at 7 < 7 ° might be analogous to such endpoints. 


7 The spectrum in the Coulomb phase 


7 .1 Expected properties of photon and resonance states 


In order to prepare the ground for an understanding of our results, let us first 
discuss the properties of a phase with massless photons and massive GB’s on 
a hnite lattice. We will conclude below that this is just what we observe from 
our data on the GB operator correlation functions. 

A massless state projected to zero (spatial) momentum does not decay ex¬ 
ponentially. From experience with Goldstone bosons [43] we expect in fact a 
propagator polynomial (parabolic) in the time variable. However, the operator 
projected to non-zero momentum will have non-zero energy E = \p\ oc 1/Ls 
and an L^-dependent exponential decay. This explicit hnite size dependence 
is an excellent indicator for massless states. Since states with non-zero mo¬ 
mentum are not parity-eigenstates, we expect that in this case the photon 
continuum state with = 1 contributes to the (1) channel. There¬ 
fore, in the ( 0 ) channel there is no signal from the photon, while the p = 1 
state has just the energy corresponding to one unit of momentum. Such obser¬ 
vations thus demonstrate the presence of a massless photon with continuum 
quantum numbers 1 . 

If, in addition, there is another, now massive, state e.g. in the Ti~~ channel, it 
will couple to multi-photon states with these quantum numbers and therefore 
will not be an asymptotic state of the system, but a resonance. The lowest 
energy state for p = 0 is the 3-photon state with photon 3-momenta (1,0,0), 
(0,1,0) and (-1,-1,0) and total energy E = {2 + \/2) (again in units of 2tt/L g). 

For simplicity of our presentation we use here the continuum dispersion rela¬ 
tion, although actually on the lattice (for massless states) we have cosh E = 
4 — cos Pi- In the actual analysis and the plots we always use the lattice 

dispersion relation. 

In a hnite volume system at hxed value of the coupling we therefore expect 
a specihc behaviour of the (discrete) energy spectrum [44] when Eg increases. 
For small Eg the lowest measured energy will be constant with a value close to 
the resonance mass and only a weak dependence on Eg. At some Eg the lowest 
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multi-photon state with energy oc l/L^ becomes the lowest energy state. The 
phenomenon of avoided level crossing will be observed, and the lowest energy 
will drop inversely proportional to the spatial extension. If one can observe 
higher levels in the energy spectrum one will see the resonance as the next 
higher state, until further avoided level crossings with multi-photon states of 
higher momenta occur. From the energy levels one can (under certain condi¬ 
tions) derive values of the phase shift in the multi-particle channels [44,33,45] 

Consider now the variation with the coupling constant. The resonance mass 
in lattice units then varies with the scale, whereas e.g. the zero-momentum 
3-photon energy in lattice units is still (2 -|- \/2). The avoided level crossing 
therefore occurs at a different lattice size. This behaviour provides another 
strong clue for the interpretation of the spectrum. 

Finally, for a given lattice size, but as a function of the coupling, the lowest 
energy level in a resonance channel might increase when moving away from the 
phase transition; eventually it should approach but never exceed the lowest 
multi-particle energy level. This level decreases proportional to l/Lg. 


7.2 Evidence for resonances 


Let us now discuss our results for gauge-balls in the Coulomb phase. Indeed, we 
find a massless state in the p = (1, 0, 0) channel of Ti'~, with a size dependence 
following the expected dispersion relation. We identify this signal with the 
photon. In the zero momentum channel we find a lowest energy state following 
the scenario of a resonance coupling to a 3-photon state, as discussed above. 
For small lattice sizes the lowest energy level is compatible with a constant, 
until it reaches the 3-photon level, which it follows for larger spatial volumes 
as seen in fig. 8. 

We estimate the position of the resonance from the energy values determined 
at lattice sizes below the point where the avoided level crossing occurs. We 
indeed find a consistent behaviour indicating an increase of the correlation 
length towards the phase transition. 

The energy levels in the channel exhibit a similar behaviour (fig. 9). This 
state is also supposedly a resonance, since it couples to a two-photon state 
(photon three-momenta (1,0,0) and (-1,0,0)) with total energy 2 (in units of 
2ttILs). 

In order to provide a further confirmation for the nature of this state we 
have also determined the second-lowest energy level at /? = 1.17 and various 
spatial volumes. As discussed above, we have usually determined the optimal 
operator from the diagonalization (cf. (17) of the correlation matrix at time 
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Fig. 8. (0) energy vs. Lg for [3 = 1.18,1.17,1.161 in the Coulomb phase. The full 

line is the energy of the 3-photon state. The dashed horizontal lines indicate the re¬ 
gion of those points that have been included in the mean value for the determination 
of the resonance positions. 

distance to + 1 = 1. Performing the diagonalization at a larger time might 
improve results for the excited states. We have tried this, but see no good 
indication that the effective mass plateau for the higher states is improved, 
and consequently have used to = 0 again. Fig. 10, albeit with larger errors and 
uncertainties for the higher lying levels, seems to support our interpretation. 

Again we determine the approximate resonance position and we hnd indication 
of critical behaviour towards the phase transition. In £g. 11 we give the j3- 
dependence for the masses of and Ti~ states. 

The data is not of sufficient quality to decide on the scaling parameters. The 
determination of masses in the Coulomb phase has various handicaps. On one 
hand, according to our observations we expect that these states are not asymp¬ 
totic but resonances. Only a narrow window is used to determine their masses. 
On the other hand, [25] observed certain effects on the mass measurements 
of the photon propagator due to the Dirac sheet background and related to 
gauge fixing [23,24]. We do not expect these to affect our results noticeably 
since we have comparatively large spatial lattices. Finally, due to the massless 
state the finite size effects are definitely larger than in the confinement phase. 



/?=1.18 


34 








2.0 


+ 




Ls 


Fig. 9. yl]^'''(0) energy vs. Lg for (3 = 1.18,1.17,1.161 in the Coulomb phase. The 
full line is the 2-photon energy. The dashed horizontal lines indicate the region of 
those points that have been included in the mean value for the determination of the 
resonance positions. 

We think that this last problem may be the main reason for our difficultiess 
to identify a consistent scaling form from the Coulomb side. 

Our results indicate that in the Coulomb phase one has both massless vector 
states (photons) and massive GB resonances and Ti~ that couple to 2 or 
3 of the massless states, respectively. The masses of the resonance states seem 
to scale towards the phase transition. We cannot decide whether both scale 
differently or not. However, assuming that the uncertainties associated with 
the determination of each individual mass partly cancel in their ratio, one can 
look at whether the log-log plot, analogous to £g. 4 in the confinement phase, 
has a slope different from one. In fig. 12 we see that the data indeed indicate 
the slope being consistent with that in fig. 4. This could point towards the 
existence of two mass scales in the Coulomb phase, too. 

As for the other channels, our data suggest that the T2*““(0), A^“(0) and 
(0) behave similar to T]^“(0). Thus massive three-photon resonances, pos¬ 
sibly with the same or similar mass, could be present in all these channels. 
The spectrum of the resonances in the Coulomb phase might thus resemble 
that of the gauge-balls in the confinement phase. 


35 
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Ls 

Fig. 10. First two energy levels of the ^^"’"(0) energy vs. Lg for beta=1.17. The 
horizontal line is the fit to the mass on small lattices. The other line represents the 
lowest two-photon energy. 

8 Conclusions and open questions 

Our results in the confinement phase are in agreement with the earlier finite 
size studies on spherical lattices [8-10] but substantially extend that work. 
They strongly suggest the existence of a continuum quantum field theory in 
four dimensions with the following properties: 

(1) From the confining phase one approaches a non-Gaussian fixed point, so 
the theory is presumably interacting but not asymptotically free. The 
correlation length exponent at this hxed point is z/ng ~ 0.36. This value 
is obtained by combining our present results with those in [8-10]. 

(2) The well measurable physical states at this hxed point contribute to 

the A 2 ,E,Ti and T 2 representations of the lattice symmetry group with 
PC = H—. They have equal or very similar mass Assuming that 

the smallest spin dominates, the continuum quantum numbers of gauge- 
balls would be = l’''“,2^ ,3^ . This assignment is uncertain as it 
could be that for example the 3’*'“ GB dominates the T^~ channel and 
the I"*"" state is not present. 

(3) Depending on the scenario for the continuum limit of the lattice theory. 
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2.0 



Fig. 11. Resonance masses vs. beta in the Coulomb phase for (0) (squares) and 
(circles). For comparison we show also masses in the same channels in the 
confinement phase. 

the continuum theory might contain a light (compared to m^g^®), possibly 
massless gauge-ball. However, it is probably only weakly cou¬ 

pled to the states in other channels, and might be completely decoupled. 
No other state in the channel was observed. 

(4) The theory is conhning in the sense of the Wilson criterium, i.e. the string 

tension in physical units has a hnite nonvanishing value — 

0.34mPjyh 

(5) In the and (presumably channels a state with 171 ^^+ ~ 
2mPg^® is present. It may be a two-particle state. 

(6) In many other channels observable states with masses ~ (2.6 — 3.6)mPg^® 
are present. Some or all of them might be two- or three-particle states. 

(7) No effective energy plateau was found in the channels Ai~,A 2 ^,E~^ 
and Ti^. However, even in these channels the effective energy at the 
distance f = 0 is measurable, indicating the possible presence of some 
heavy, possibly multiparticle states. 

We point out that because of the rigorous duality relations for the U(l) lat¬ 
tice gauge theory the same continuum theory can also be obtained as the 
continuum limit of the following four-dimensional lattice theories: Coulomb 
gas of monopole loops [46], Z gauge theory [47], and non-compact U(l) Higgs 
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Fig. 12. Different scaling of the two gauge-ball states observed in the Coulomb 
phase, assuming that the systematic uncertainties are less important in the ratio of 
the two masses. 


model at large negative squared bare scalar mass (frozen four-dimensional su¬ 
perconductor) [2,47]. The last model can be understood as a limit of a theory 
described by world sheets of Nielsen-Olesen strings [48]. 

Currently we have no physical picture explaining the observed gauge-ball spec¬ 
trum. It could be that its understanding might rather come from one of the 
dual equivalents, for example as the lowest states of a closed Nielsen-Olesen 
string. 

Our results in the Coulomb phase indicate that also here a continuum theory 
approaching the conhnement transition might exist which differs from the free 
Maxwell theory usually expected in this phase at least in the /? —> cx) limit. 
However, the masses and the scaling behaviour of the resonances observed in 
this phase are not yet fully understood. One or even two mass scales might 
emerge there. 

As a spin-off of our study of this phase we conhrm the universality [13-15] of 
the maximal value of the renormalized coupling by obtaining = 0.19(1) in 
agreement with results in other lattice formulations of the U(l) gauge theory. 
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Our present work leaves unanswered several questions about the pure U(l) 
lattice gauge theory and raises many new ones. The major challenge is to un¬ 
derstand how the order of the conhnement-Coulomb phase transition changes 
from hrst at 7 > 70 to second at 7 < 79. Is the point 70 a tricritical point in 
the sense that several critical lines emerge from it? What are the additional 
couplings required to reveal this structure? Where are the crossover regions 
between critical and tricritical scaling behaviour? Are these regions for the O’*"’*' 
state different from those of the other GB states, as required by the Tng and 
Tg scenarios? Or is 70 an ordinary point of the critical line 7 < 7o, as would 
be most natural in the scenario C? Is z/^g a critical or tricritical exponent? A 
clarification of these questions will require large effort, as 7 and possibly some 
other couplings will have to be varied. 

Some arguments can be given for prefering the scenario Tng: (i) It is consistent 
with some MCRG studies [12]. (ii) Previous results for the string tension [15] 
and for the massive photon mass [23] at 7 = 0, i.e. very close to 70, are 
consistent with the non-Gaussian value, (hi) It is a fairly conventional scenario 
implying a broad range of dominance of the TGP. 

However, the observation of z/ ~ 0.36 in the hnite size scaling analysis even 
at 7 = —0.5 [8,9] could be difficult to accommodate in that scenario and 
prefers scenario C. Two different scaling laws at the same critical point might 
seem to be somewhat exotic. However, the rigorously established results in 
the three-dimensional pure U(l) gauge theory [37] are strikingly similar. Also 
there a 0 ’*'^ mass scales faster than the string tension when the gauge coupling 
vanishes (both scale exponentially). This implies the existence of two different 
continuum limits. The one with conhnement contains a massless scalar, too. 

We stress that the existence of the non-Gaussian continuum theory does not 
depend on which scenario is correct. Only the properties of the state, 
whether it is light but massive or massless, and whether it couples or not to 
the other states, depend on the scenario. Thus the fate of the state in the 
non-Gaussian continuum limit is a challenging question. 

The particular scaling behaviour of the GB observed in the present work, 
in the lattice QGD [38] and in the three-dimensional pure U(l) lattice gauge 
theory [37] should draw attention towards light scalars in the conhning theo¬ 
ries. 

In the light of so many open questions we should mention the natural caveat 
of numerical simulations: it might be that only much closer to the phase tran¬ 
sition, on lattices much larger than we were able to use, the true behaviour 
will show up in the future. For example, the scaling of all observables might 
turn out to be the same at any 7 < 70, or the first order transition might 
reappear even on spherical lattices. 
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There are several paths which could be pursued with the means currently 
available: 

• Studies similar to the present one at various 7 could elucidate the fate of 
the state and show the cross-over regions if they exist. For this purpose 
both large negative and some positive 7 values will be required. 

• It should be possible to investigate the short-range form of the static poten¬ 
tial in the conhnement phase. This would shed more light on the nontrivial 
held theory strongly interacting at short distances. 

• It would be interesting to hnd out whether the heavier states observed in 
the conhnement phase are genuine gauge-balls or multi-particle states. If 
they are resonances, they should be investigated with the appropriate hnite 
size techniques [44,33,45]. 

• An analysis of the resonances in the Coulomb phase with these techniques 
would establish their existence and allow to determine the scaling behaviour 
and continuum limit in this phase as well. The obstacle is the presence of the 
massless photon, which complicates the analysis devised for theories with a 
mass gap [44]. 

• Simulations of the dual Z theory might bring further insight into the £- 
nite size properties, scaling behaviour of the monopole mass [26], role of 
boundary conditions, etc. 

The pure gauge theory investigated here may not be realized in nature. How¬ 
ever, it seems worthwhile to pursue its study, since it widens our understanding 
of quantum field theories. Our lattice simulations have brought about unex¬ 
pected results. Further surprises may be waiting. 
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A Appendix 


In this appendix we present two tables with more detailed information about 
the GB masses. Further data on the effective energies in the conhnement phase 
can be obtained from the authors by e-mail (jersak@physik.rwth-aachen.de). 

Table A.l 

The effective energies of the and gauge-balls with zero momentum in the 
confinement phase on the time distances t/{t + 1). 


A++ 

Gitter 

0 

0/1 

1/2 

2/3 

3/4 

4/5 

4/6 

6/7 


1.13 

1.159(8) 

1.09(2) 

1.05(6) 

0.90(16) 

1.1(7) 

1.3(11) 



1.135 

1.074(8) 

1.023(20) 

0.96(7) 

1.05(15) 

1.2(8) 

1.4(11) 



1.14 

0.968(6) 

0.914(14) 

0.92(3) 

1.05(11) 

0.7(3) 

0.4(3) 

0.2(5) 


1.145 

0.862(6) 

0.815(14) 

0.84(4) 

0.88(8) 

0.76(16) 

0.7(3) 

0.1(2) 

16^32 

1.15 

0.732(4) 

0.692(9) 

0.701(15) 

0.73(4) 

0.77(9) 

0.7(2) 

0.4(4) 


1.152 

0.654(6) 

0.618(10) 

0.595(15) 

0.58(3) 

0.60(5) 

0.60(11) 

0.41(13) 


1.154 

0.607(8) 

0.559(12) 

0.557(18) 

0.53(3) 

0.50(5) 

0.47(9) 

0.42(15) 


1.156 

0.511(7) 

0.472(10) 

0.468(15) 

0.46(2) 

0.48(3) 

0.46(5) 

0.50(8) 


1.158 

0.412(10) 

0.375(12) 

0.376(16) 

0.38(2) 

0.39(3) 

0.45(5) 

0.48(8) 


1.154 

0.588(10) 

0.544(13) 

0.530(18) 

0.53(3) 

0.57(6) 

0.68(17) 

1.1(8) 


1.156 

0.533(11) 

0.490(17) 

0.50(3) 

0.51(5) 

0.54(9) 

0.54(19) 

0.5(3) 

20^40 

1.158 

0.410(10) 

0.365(14) 

0.355(19) 

0.35(3) 

0.33(4) 

0.30(5) 

0.28(6) 


1.159 

0.306(5) 

0.271(6) 

0.268(7) 

0.262(8) 

0.263(10) 

0.262(12) 

0.260(16) 


T+- 

Gitter 

0 

0/1 

1/2 

2/3 

3/4 

4/5 

4/6 

6/7 


1.13 

1.496(6) 

1.43(2) 

1.32(12) 

1.5(4) 





1.135 

1.384(5) 

1.314(20) 

1.33(8) 

1.1(3) 

1.4(11) 




1.14 

1.293(5) 

1.236(16) 

1.24(6) 

1.17(17) 

0.6(4) 

0.3(5) 

1.4(9) 


1.145 

1.189(4) 

1.121(12) 

1.08(3) 

0.99(10) 

0.9(2) 

0.4(4) 

0.4(5) 

16332 

1.15 

1.066(4) 

0.981(9) 

0.97(3) 

0.92(6) 

0.80(15) 

0.6(3) 



1.152 

1.006(4) 

0.925(8) 

0.903(19) 

0.90(6) 

0.94(16) 

1.1(9) 



1.154 

0.950(4) 

0.863(8) 

0.843(18) 

0.87(4) 

0.87(12) 

1.2(5) 

0.4(7) 


1.156 

0.798(4) 

0.738(8) 

0.732(14) 

0.75(3) 

0.73(5) 

0.71(14) 

1.0(5) 


1.158 

0.668(4) 

0.601(6) 

0.588(10) 

0.564(17) 

0.56(3) 

0.57(6) 

0.54(10) 


1.154 

0.956(5) 

0.862(13) 

0.83(3) 

0.76(5) 

0.73(15) 

0.5(2) 

0.6(6) 


1.156 

0.779(6) 

0.729(11) 

0.74(3) 

0.73(5) 

0.81(11) 

0.7(2) 

0.8(5) 

20^40 

1.158 

0.676(5) 

0.624(9) 

0.624(18) 

0.60(3) 

0.60(5) 

0.58(11) 

0.7(2) 


1.159 

0.5933(19) 

0.528(3) 

0.518(4) 

0.510(6) 

0.528(11) 

0.540(16) 

0.53(3) 
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Table A.2 

Masses for all combinations obtained from operators with zero and smallest 
nonvanishing momentum. 


RPC 

VzL 

27r 

16^32 

/3 = 1.13 

/3 = 1.135 

/3 = 1.14 

/3 = 1.145 

/? = 1.15 

/3 = 1.152 

P = 1.154 

A++ 

0 

1.08(7) 

1.02(8) 

0.95(4) 

0.85(5) 

0.721(20) 

0.599(19) 

0.55(2) 

1 

1.27(11) 

1.25(8) 

1.00(4) 

0.80(3) 

0.70(3) 

0.66(2) 

0.566(20) 


0 



2.3(12) 

3.9(18) 

2.1(10) 

1.8(4) 

2.1(6) 


0 

2.6(17) 




2.4(6) 

2.7(17) 

2.7(8) 


0 

3(2) 


1.0(11) 

1.4(12) 

2.1(7) 

1.8(14) 

2.1(7) 


0 

2.2(14) 

2.5(16) 


5(4) 

3.5(19) 

5(3) 

1.9(4) 

At~ 

0 

1.5(3) 

1.6(3) 

1.20(10) 

1.19(9) 

1.00(4) 

0.95(5) 

0.92(4) 


0 

1.3(13) 



1.9(11) 

3(2) 

5(3) 


^2 

0 


3.2(20) 


11(8) 

3.6(17) 


1.2(6) 

E++ 

0 

2.2(5) 

1.8(9) 

3.0(5) 

2.6(3) 

2.10(14) 

1.88(12) 

1.74(8) 

1 

1.4(14) 

2.8(7) 

2.5(3) 

2.6(3) 

2.25(17) 

2.04(13) 

1.99(11) 

E+- 

0 

1.65(19) 

1.29(12) 

1.31(11) 

1.05(6) 

0.96(4) 

0.91(3) 

0.84(3) 

E-+ 

0 


6(2) 

1.2(16) 

1.9(9) 

4.0(16) 

2.9(6) 

3.0(7) 

E — 

0 

10(4) 

7(5) 


5(3) 

2.2(6) 

2.3(6) 

2.7(6) 


0 


4(2) 

10(5) 


4(2) 

9(4) 

1.4(9) 

1 


5(2) 



2.3(12) 


1.4(5) 

T+- 

0 

1.43(15) 

1.42(9) 

1.29(7) 

1.11(4) 

1.00(3) 

0.93(2) 

0.87(2) 

1 

1.58(19) 

1.25(16) 

1.50(12) 

1.16(9) 

1.01(4) 

0.85(4) 

0.88(4) 

T-+ 

0 

2.4(15) 

2.5(15) 

1.4(8) 

5(3) 

13(5) 

1.8(9) 

4(2) 

1 

1.7(4) 

1.0(3) 

1.13(11) 

0.76(7) 

0.73(3) 

0.67(4) 

0.59(3) 

T" 

0 

2.6(17) 

4(2) 

30(13) 

1.6(8) 

2.7(11) 

1.7(7) 

3.9(15) 

1 

0.8(14) 

1.4(6) 

1.6(11) 


2.7(14) 

4(2) 

1.7(14) 

^2 

0 

3.6(11) 

3.3(6) 

3.1(3) 

2.6(2) 

2.13(11) 

1.91(8) 

1.85(6) 

1 

4.8(20) 

2.5(4) 

4.4(16) 

2.3(2) 

2.2(2) 

1.95(14) 

1.73(8) 

jiH — 
^2 

0 

1.48(14) 

1.44(10) 

1.25(8) 

1.07(5) 

0.99(3) 

0.93(2) 

0.86(2) 

1 

1.37(19) 

1.14(14) 

1.45(14) 

1.28(9) 

0.99(5) 

0.92(5) 

0.83(4) 

rp -h 

^2 

0 

6(3) 

2.6(12) 

2.3(10) 

1.7(14) 

3.9(12) 

2.4(4) 

2.1(2) 

1 



1.7(4) 

1.8(13) 

3.0(9) 

2.4(6) 

2.4(4) 

T2~ 

0 



7(3) 

19(10) 

7(4) 

2.8(17) 

2.8(7) 

1 

2.9(19) 

3(2) 

1.0(4) 

1.4(7) 

0.9(2) 

1.09(15) 

0.86(9) 
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Table A.3 

Continuation of table A.2 


RPC 

PzL 

27r 

10332 

20^40 

l3 = 1.156 

/? = 1.158 

/? = 1.154 

/? = 1.156 

f3 = 1.158 

l3 = 1.159 

A++ 

0 

0.47(2) 

0.39(2) 

0.54(2) 

0.51(4) 

0.34(3) 

0.265(9) 

1 

0.459(19) 

0.33(3) 

0.51(2) 

0.43(3) 

0.31(2) 

0.293(6) 


0 

1.88(18) 

1.70(11) 

2.3(7) 

1.8(9) 

1.6(2) 

1.65(4) 


0 

2.7(6) 

1.72(14) 

11(6) 

1.9(6) 

2.0(3) 

1.49(3) 


0 

2.6(7) 

1.73(20) 

1.6(14) 

1.2(3) 

1.3(4) 

1.47(4) 


0 

2.3(13) 

2.0(2) 

4(2) 

2.2(8) 

2.1(3) 

1.70(5) 

At~ 

0 

0.79(3) 

0.61(2) 

0.84(6) 

0.71(4) 

0.62(3) 

0.529(10) 


0 


1.5(4) 

1.0(17) 

0.8(9) 


2.8(4) 

^2 

0 

2.4(9) 

1.8(2) 

1.9(12) 

2.1(10) 

1.8(4) 

1.46(5) 

E++ 

0 

1.56(5) 

1.22(2) 

1.70(9) 

1.60(8) 

1.30(5) 

1.098(11) 

1 

1.67(7) 

1.47(4) 

1.90(12) 

1.80(13) 

1.42(6) 

1.193(12) 

E+- 

0 

0.76(2) 

0.590(17) 

0.85(4) 

0.76(4) 

0.62(3) 

0.510(8) 

E-+ 

0 

2.1(3) 

1.63(10) 

1.9(13) 

1.5(2) 

1.66(15) 

1.58(3) 

E — 

0 

2.02(19) 

1.78(10) 

1.6(4) 

2.1(4) 

1.77(15) 

1.56(3) 

T++ 

0 

0.9(7) 

2.4(3) 

2.8(15) 

10(5) 

1.8(14) 

1.94(6) 

1 

3.0(11) 

1.9(3) 

6(3) 

1.4(8) 

2.9(12) 

1.89(9) 

T+- 

0 

0.751(16) 

0.588(12) 

0.83(3) 

0.76(4) 

0.63(2) 

0.524(5) 

1 

0.76(3) 

0.64(2) 

0.78(5) 

0.71(3) 

0.62(2) 

0.522(6) 

T- + 

0 

1.8(2) 

1.9(4) 

2.5(14) 


1.6(3) 

2.22(9) 

1 

0.45(2) 

0.33(3) 

0.51(3) 

0.45(4) 

0.33(2) 

0.295(8) 

T” 

0 

1.81(17) 

1.62(9) 

3.0(15) 

1.9(4) 

1.82(16) 

1.50(3) 

1 

1.7(2) 

2.09(18) 

2.3(8) 

1.6(8) 

1.7(2) 

1.59(3) 

rp + -\- 
^2 

0 

1.68(5) 

1.29(3) 

1.75(9) 

1.53(6) 

1.26(4) 

1.077(8) 

1 

1.64(6) 

1.34(4) 

1.66(12) 

1.55(10) 

1.31(5) 

1.209(12) 

rp-\-- 

0 

0.761(16) 

0.591(12) 

0.85(3) 

0.76(4) 

0.63(2) 

0.525(7) 

1 

0.75(3) 

0.60(2) 

0.88(6) 

0.75(5) 

0.61(3) 

0.533(9) 

rp -h 

-^2 

0 

1.68(13) 

1.78(8) 

2.2(3) 

1.9(3) 

1.79(13) 

1.69(3) 

1 

1.87(20) 

1.59(9) 

1.8(3) 

1.8(4) 

1.42(11) 

1.38(2) 

rp-- 

-^2 

0 

1.93(20) 

1.99(14) 

4(2) 

4(2) 

1.91(19) 

1.74(5) 

1 

0.84(6) 

0.65(4) 

0.68(16) 

0.63(7) 

0.71(6) 

0.548(13) 
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